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Abstract 

Some recent progress in the mathematical physics of rapidly rotating, dilute Bose gases in anhar- 
monic traps is reviewed. 

1 Introduction 

Since the first experimental realization of Bose-Einstein condensation in trapped, dilute alkali gases in 
1995 the strange quantum properties exhibited by ultracold quantum gases have become a large research 
field worldwide, both experimentally and theoretically. Under normal conditions the atoms of a dilute 
gas are distributed on very many one-particle quantum states so that every single state is only occupied 
by relatively few atoms. By cooling a gas of alkali atoms to extremely low temperatures (of the order 
10~ 8 K) it is, however, possible to achieve Bose-Einstein condensation, where a single quantum state 
is, on the average, occupied by a macroscopic number of particles. This phenomenon was predicted by 
Albert Einstein in 1924 based on a model of an ideal quantum gas, where the unavoidable interactions 
between the atoms are ignored. The interactions, however, are important and modify the phenomenon in 
several ways. Despite much progress in the theoretical understanding (see, e.g., the monographs pQ and 
[2]) it is still a real challenge for mathematical physics to derive the properties of the low energy states 
of interacting Bose gases from the many-body Hamiltonian by rigorous mathematical analysis. I refer to 
[3] for an introduction to these problems. 

In the present contribution the focus will be the effects of rotations on the ground state properties 
of dilute Bose gases. In particular, I shall review some recent results obtained in collaborations with 
Jean-Bernard Bru, Michele Correggi, Tanja Rindler-Daller and Peter Pickl [H El El I7J • These results 
concern mainly the ground state energy of a gas in a rapidly rotating container as well as the vortices 
that are created by the rotation and arrange themselves in a geometric pattern to minimize the energy. 
The external potential modelling the container is here assumed to increase faster than quadratically with 
the distance from the rotational axis ('anharmonic traps'), allowing arbitrarily high rotational velocities. 
By contrast, in a quadratic confining potential ('harmonic trap') the centrifugal forces will drive the gas 
out of the container when the rotational velocity exceeds a finite value. When this limiting velocity is 
approached from below interesting new effects are expected to occur but they are not subject of the 
present review except for some brief comments in Section [7] 
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2 The General Setting 

We consider the quantum mechanical Hamiltonian for N bosons with a pair interaction potential, v, and 
external potential, V, in a rotating frame with angular velocity ft: 

N 

F = ]T(-V J 2 + y(x J )-L J .fi)+ Yl (2- 1 ) 

j=l l<i<j<N 

Here x^- G R 3 , j = 1, . . . , N are the positions and Lj = — ix^ A Vj the angular momentum operators of 
the particles. Units have been chosen so that ft = 2m = 1. 

The pair interaction potential v is assumed to be radially symmetric and of finite range. For the 
mathematical results in the sequel we also assume that v is nonnegative. The Hamiltonian operates on 
symmetric (bosonic) wave functions in L 2 (R 3N , d 3 xi • • • d 3 XAr). 

The Hamiltonian can alternatively be written in the form 

N 

H = J2((-™ ] -Mxj)f + V(x 1 )-\tfr])+ ]T v(\ Xi - Xj \) (2.2) 

j = l l<i<j<N 

with A(x) = ifiAx and r the distance from the rotation axis. This corresponds to a splitting of the 
rotational effects into Coriolis and centrifugal forces. The Coriolis forces have formally the same effect 
as a magnetic field with vector potential A(x). 

We are interested in the ground state properties of the Hamiltonian for N — > oo, in particular the 
effect of rotations on dilute gases. Dilute means that 

a 3 ( o<l (2.3) 

where a is the scattering length of v and p an appropriately defined average particle density. The scattering 
length is defined by considering the zero energy scattering equation 

- V 2 ?A + \vil> = 0. (2.4) 

For |x| larger than the range of v the solution is 

V>(x) = (const.) (l - alxl" 1 ) (2.5) 

which defines a. The condition a 3 p <C 1 can also be written 

a«p- 1/3 (2.6) 

i.e., the scattering length is much smaller than the mean particle distance. 

A special case of a dilute limit is the Gross-Pitaevskii (GP) limit, which means that the GP coupling 
parameter 

g = AiraN/L (2.7) 

is kept constant as N — > oo. Here L is the length scale associated with —V 2 + V, measured, e.g., in terms 
of the width Ax of the ground state of this operator. Note that p <~ N/L 3 , so in the GP limit 

a 3 p~g 3 /N 2 = 0(N- 2 ). (2.8) 

To understand the meaning of the GP limit one should recall a basic formula for the ground state 
energy per particle of a dilute Bose gas in a box: 

E a /N w Airpa. (2.9) 
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Thus, in the GP limit, the interaction energy per particle, ~ Na/L 3 , is of the same order as the lowest 
excitation energy without interaction, ~ l/L 2 . 

The formula (|2.9[) has an interesting history, spanning almost 80 years! The technique used for the 
first complete mathematical proof [5] (see also [H] for the two-dimensional case) has been important for 
subsequent rigorous work on dilute Bose gases. 

Associated with the energy ~ pa there is a natural length scale, the 'healing length' 

4 = (Anpa)- 1/2 . (2.10) 

Thus 1/i 2 . is equal to the interaction energy per particle. If the density of an inhomogeneous gas is 
somewhere zero, due to boundary conditions or the formation of vortices, the interaction will tend to 
restore the full density ('heal') on a length scale of order £ c . Note that £ c decreases if the scattering 
length increases. For later use it is convenient to define 

e = e c /L = g- 1 / 2 . (2.11) 

In rotating gases e is a measure for the size of the vortex cores relative to the size of the trap. 
Besides the energy it is also of interest to consider the particle density 

po(x)=Jv/ |* (x,x 2 ,...,x JV )| 2 d 3 x 2 ---d 3 x A r (2.12) 

associated with a (ground state) wave function ^!q, and the one-particle density matrix 

7o (x,x') = W *o(x,x 2 ,...,x Ar )^(x',x 2 ,...,x A ,)d 3 x 2 ---d 3 x JV . (2.13) 

Bose-Einstein condensation (BEC) in the state means (by definition) that the operator defined by 
the integral kernel 70 (x, x') has an eigenvalue of the order N for all large N. This means that 

7o (x,x') = A Vo(x)MxT +Rest (2.14) 

with Ao = O(N) for large N. The ratio Xq/N is called the condensate fraction and tpa the wave function 
of the condensate (or sometimes the 'superfluid order parameter'). The 'rest' in (|2.14[) corresponds to 
the depletion of the condensate, i.e., the contribution of states orthogonal to ho- 
using creation and annihilation operators, we can more generally write 

7 o(x,x')-(a t (x) a (x')> (2.15) 

where (•) is the expectation value in the many-body state that need not be a pure state and could, e.g., 
be a thermal equilibrium state. Thus 

Ao = (a^oM^o)) (2.16) 

is the average particle number in the one-particle state tpo when the whole system is in the given many- 
body state. 

The following basic results on the GP limit of the ground state of the many-body Hamiltonian for 
g = AnNa/L and f2 fixed were established in [TU 1 fTT | fT2"]: 

• The ground states can be described through minimizers of the Gross- Pitaevskii energy functional 

£ GP [V>] = / {|v^| a + ^|Vl a -n.v*i^+ffH 4 }d 3 x 

JR 3 

= / {|(iV+ A)^| 2 + (y-ifiV)M 2 +. 9 h/f }d 3 x (2.17) 

JR3 

with the normalization condition /rsI^I 2 = 1. 
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• For rotating gases the bosonic ground state is in general not the same as the absolute ground state 
and it need not be unique due to breaking of rotational symmetry. 

More precisely, Bose-Einstein condensation takes place in the GP limit and every one particle density 
matrix, obtained as the limit of one-particle density matrices, divided by N, of ground states of (|2.2p is 
a convex combination of projectors onto minimizers of the GP energy functional. 

The GP minimizers are solutions of the GP equation 

{-(V - i A) 2 + (V — |0 2 r 2 ) + 2g\ib\ 2 } V = (2.18) 

It has the form of a non-linear Schrodinger equation, with a magnetic field in the case of rotation. The 
main new feature compared to the non-rotating case is the possible occurrence of vortices, i.e., singularities 
of the phase of ip with integer winding numbers. 

The GP equation and its vortex solutions is a subject of its own that can be studied independently 
of the many-body problem. See the monograph [T3] and the review article [H] where a large number of 
references can be found. The most detailed results are for the two-dimensional GP equation, i.e., when ip 
depends only on the coordinates in the plane perpendicular to the angular velocity f2. In particular, the 
two-dimensional GP equation with a quadratic trap potentials V^(r) ~ r 2 has been studied in P^lfTBlfTT] 
in the asymptotic parameter regime 

g^co, Q ~ g-l/ 2 \l gg\. (2.19) 
More general, homogeneous trapping potentials are discussed, e.g., in [51 118] . 



3 Two- Dimensional GP Vortices 

This section provides some heuristic background for understanding the occurrence of vortices in the case 
of the two-dimensional GP equation. In order to bring out the salient points as simply as possible we 
consider a 'flat' trap, i.e., the condensate is confined to a circle of radius R and we impose Neumann 
conditions at the boundary. 

The first thing to note is that for sufficiently small rotational velocities the condensate stays at rest in 
the inertial frame and thus appears to rotate opposite to CI in the rotating frame. This is a manifestation 
of superfluidity: A normal fluid would pick up the rotational velocity of the container and in equilibrium 
the fluid would be at rest in the rotating frame. 

In the rotating frame the operator of the velocity is — iV — A(r). The constant wave function, that 
minimizes the GP energy functional (with zero energy in excess of the interaction energy) for small f2, 
thus has, in the rotating frame, the velocity 

v(r) = -A(r) = -|0 Ar = -§Qre e , (3.1) 

where eg denotes the unit vector with respect to the angular variable. Note that the kinetic energy 
corresponding to this velocity is exactly compensated by the centrifugal term — if2 2 r 2 in the GP energy 
functional pT7]) , 

At higher rotational velocities the condensate responds by creating vortices in order to reduce the 
velocity. Consider the case of large g, i.e., small e ~ g~ x ^ 2 ■ A vortex of degree d located at the origin has 
the form 

^(r,0) = /(r)exp(i0d) (3.2) 

with 

( r d if < r < eR 
fir) ~ { (3.3) 
R- 1 if sR<r<R 
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Now the component of the velocity in the direction of eg is 

v(r) e = (i - ISlA eg . (3.4) 
The energy (in excess of the interaction energy) is therefore 

~ R- 2 [{d/rf -dn]rdr + 0{\) = R- 2 d 2 \\oge\-\dtt + 0(l). (3.5) 

JeR 

Thus a vortex of degree d = 1 can be expected to form when 

ir 2 |loge| - in < (3.6) 

i.e., 

n > Oi - ir 2 |loge|. (3.7) 

We also see that d vortices of degree 1, ignoring their interaction, have energy ~ d(R~ 2 \ loge| — fi/2) 
while a vortex of degree d has energy R~ 2 d 2 \ log e| — dil/2. Hence it is energetically favorable to 'split' 
a d-vortex into d pieces of 1-vortices, breaking the rotational symmetry. 

These heuristic considerations are confirmed by a detailed analysis for 'slowly' rotating gases, i.e., 
fl = 0(| logeQ [mUSHHHH]: Vortices start to appear for ttR 2 = n\ loge\ and for 

|loge| + (d-l)log|loge| < flR 2 /n < | loge| + d\og | loge| (3.8) 

there are exactly d vortices of degree 1. If the 'flat' trap is replaced by a homogeneous trap, V(r) ~ r s , the 
effective radius R of the condensate is determined by equating the potential energy and the interaction 
energy, i.e., R s ~ e^ 2 R 2 (R^ 1 ) 4 , which leads to 

R ~ £ - 2 '^ (3.9) 

and thus (|3.8p is replaced by 

Qx ~ £ 4 /< s+2 ^|loge|. (3.10) 

In |19j a number of general results about two-dimensional GP vortices are derived that are not 
limited to the asymptotic regime e — > 0. The trap potential is assumed to be rotationally symmetric and 
polynomially bounded; for V(r) ~ r 2 one must in addition assume that is less than a critical velocity, 
fi c , so that V(r) — fl 2 r 2 /4 is bounded below. Among the results proved in [TH] are 

• Instability for large d: For all < £1 < f2 c there is a dn, independent of the coupling constant g, so 
that all vortices with d > do, are unstable. 

• Symmetry breaking: For all < O < O c there is a <?q, so that g > gn implies that no ground state 
of the GP energy functional is an eigenfunction of angular momentum. 

The limit g — > oo, i.e., e — > 0, is usually referred to as a Thomas-Fermi (TF) limit because the 
resulting energy functional, which is just the GP functional without the kinetic energy term, has a formal 
similarity with the one studied in TF theory for electrons although the physics is quite different. 

4 Rapid Rotation and Vortex Lattices 

This section contains a resume of [7]. The setting is again that of two-dimensional GP theory in a 'flat', 
two-dimensional trap trap and by simple scaling we may take R = 1. If fl = 0(\ loge|) there is a finite 
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number of vortices, even as e — * 0, as mentioned in the previous section. For larger f2 the number of 
vortices is unbounded as e — » 0. Of particular interest is the case O = 0(\/e). 
For the flat, 2D trap the GP energy functional is 

S GP m= [ {|(iV + A)7/f- iftV|?/f +e- 2 H 4 }d 2 r (4.1) 

J\r\<l 

with A(r) = ^ilreg. The special significance of the case ft — 0(1/ e) can be seen from the GP en- 
ergy functional: The centrifugal term — (f2 2 /4)r 2 |^>(r)| 2 and the interaction term (1/e 2 ) |-0(r) | 4 are then 
comparable in size. 

The 'kinetic energy' term |(iV + A(r))?/>(r)| 2 is formally also of order 1/e 2 if O ~ 1/e. However, 
it turns out that its contribution to the energy is, in fact, of lower order because a lattice of vortices 
emerges as e — > 0. The velocity field generated by the vortices compensates partly that generated by 
A(r) = |f2 A r. The precise statement, proved in [7], that puts the analysis of (5D] on a rigorous basis 
and generalizes previous results of [3], is as follows: 



Theorem 4.1 (Energy asymptotics, 2D GP theory) Let E GP denote the GP energy, i.e., the min- 
imum of the GP energy functional &4-l\ ). Let E TF denote the minimal energy of the GP functional without 
the kinetic term. 

If 1/e < ft< l/(e 2 |loge|) then 

E GF =E TF + ±Q\ loge| (l + o(l)). (4.2) 

// | loge| < Q, < 1/e, then 

e gp = e tf + i n j i og ( £ 2 n) | (1 + o(1)) _ (4 3) 

In particular the contribution of the kinetic term is only 0(\ loge|/e) instead of 0(l/e 2 ) if Q ~ 1/e. 
The energy E TF is the infimum of the TF functional 

£ TF [p]= I {-^Vp + e-V}d 2 r (4.4) 

J\v\<l 

where the infimum is over all nonnegative p(r) with J p = 1. The energy E TF and the (unique) minimizer, 
p TF , are easily computed explicitly. Due to the centrifugal forces the density has a parabolic shape for 
fi < Whdie/e with Whole = 4/y / 7t, but develops a 'hole' if fl > ^hoic/e- For 3> Whoio/e the density is 
concentrated in a thin layer of thickness ~ efl^ 1 at the boundary of the trap. 

In the proof of Theorem 14. II one uses, for the upper bound, a trial function corresponding to a lattice 
of vortices. Writing points r = (x,y) £ R 2 as complex numbers, £ = x + iy, the trial function has the 
form 

V>(r) = f(r)s/p~(rj exp{i<p(r)} (4.5) 
where p(r) is a suitably regularized version of the TF density p TF (r) and 

exp{V(r)} = n^-|| (4.6) 

where the points Q form a regular lattice within the trap. The function /(r) has a zero of order 1 at each 
of the points Q and is constant outside a disc ('vortex core') of radius t <C 51 -1 / 2 around each lattice 
point. 
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Note that 

p(r)=£arg(C-Ci) (4.7) 

i 

is a harmonic function. Its conjugate harmonic function is 

x(r) = ^log|r-r i |. (4.8) 

i 

Indeed, \ + iip = log(£ — Q) is an analytic function. From the Cauchy-Riemann equations 

V r x = V 9¥ > , V eX =-\7 r ip (4.9) 

it follows that 

\{V V -\nre e )\ 2 = |(V X - i^re,.)! 2 . (4.10) 

The field 

E(r) = V%- |fire r (4.11) 

can be interpreted as an electric field generated by unit 'charges' at the positions of the vortices, r^, and 
a uniform charge distribution of density f7 /2tt. 

We now distribute the vortices also with the same average density as the uniform charge. Thus 
every vortex sits at the center of lattice cell Qi of area \Qi\ = 2ir/n, surrounded by a uniform charge 
distribution of the opposite sign so that the total charge in the cell is zero. If the cells were disc-shaped 
Newton's theorem would imply that the 'electric' field generated by the cell would vanish outside the cell, 
i.e, there would be no interaction between the cells. Now the cells can not be disc shaped, but the closest 
approximation are hexagonal cells, corresponding to a triangular lattice of vortices, giving the optimal 
energy. The interaction between the cells, although not zero, is small because the cells have no dipole 
moment. 

The energy of each cell, taking into account the factor /(r) that cuts off the Coulomb field from the 
point charge (vortex) at a radius t <C fi" 1 / 2 , is 

2tt^ (1/r) 2 rdr + 0(1) = tt| log(i 2 ft)| +0(1). (4.12) 

We now multiply by the number of cells, ri/2, and divide by the area, 7r, of the trap, obtaining 
log(i 2 51)|(l + o(l)). The choice of the vortex radius t is dictated by the error terms. The opti- 
mal choice is 

• if O < 1/e 

£ -p~ 1/2 = \ (e / n) v» if i/ e < n < 4 - 13 ) 

Indeed, for f2 < 1/e the average density p is 0(1), while for £1 ^ 1/e the centrifugal forces cause a 
concentration of the GP wave function in a thin layer of thickness 0((ef2) -1 ) at the boundary and the 
average density is 0(e£l). 

A lower bound matching the upper bound to subleading order can for 

|loge| <ft< l/(e 2 |loge|) (4.14) 

be proved making use of techniques from Ginzburg-Landau theory developed by S. Serfaty and E. Sandier 
[21| I22j . An important ingredient is a splitting of the domain into small cells and a rescaling of lengths 
in each cell. This allows to reduce the problem to a study of a Ginzburg-Landau functional in small cells 
with a new interaction parameter e and a magnetic field h cx <C | loge|. 
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Outside the parameter range (|4.14[) the ansatz (|4. 5[) does not give the correct energy to subleading 
order. In fact, for f2 < | log e | there are only finitely many vortices (see, e.g., [H]-[TH]) and the upper 
bound (|4.3[) is too large. For l/(e 2 |loge|) < fi, on the other hand, a trial function different from (|4.5[) 
gives lower energy than (|4.2p . This is the trial function considered in 4 Eq. (3.36) for ^> e~ x that 
corresponds to a giant vortex [20] where all the vorticity is concentrated at the center while the support 
of p TF is vortex free. More explicitly, the trial function has the form 

^(r) = V^exp(i[n/2]0) (4.15) 

where p(r) is a regularization of the density p TF and [0/2] denotes the integer part of ft/2. With 
this trial function the next correction to the TF energy is 0(l/e 2 ) that is smaller than f2|loge| when 
l/(e 2 | loge|) < ft. This transition can be understood by the following heuristic argument, employing the 
electrostatic analogy: For CI >• 1/e the number of cells in the support of p TF is of the order 1/e. Without 
vortices each cell has unit 'charge', originating from the vector potential, and the mutual interaction 
energy of the cells is of the order 1/e 2 . Putting a vortex in each cell neutralizes the charge so that the 
interaction energy becomes negligible, but instead there is an energy cost of order fi|logej due to the 
vortices. Equating these two energies leads to ft ~ l/(e 2 | loge|) as the limiting rotational velocity above 
which the ansatz (|4.5|) is definitely not optimal. 

It should be noted that also for 1/e < ft -C l/(e 2 | loge|) one could for the upper bound replace the 
distribution of the vorticity on a lattice within the 'hole' by a single phase factor corresponding to a giant 
vortex at the origin, but in order to obtain the correction beyond the TF term the support of p TF can 
not be vortex free. The detailed vortex distribution of the true minimizer of the GP energy functional is 
still an open question. 

5 Techniques for Deriving the GP Equation 

This section contains a brief sketch of methods that have been used for a rigorous derivation of the 
three-dimensional GP energy functional from the many-body Hamiltonian. 

5.1 Nonrotating trap 

The first step is a derivation of the ground state energy formula (|2.9|) in the case of a Bose gases in a 
box. For an upper bound [S3] (see also [TUl [21]) one can use a trial function of the form 

F(xi, . . . ,xjv) = F 2 (xi,x 2 ) • ■■F N (xi, . . . ,x N ) (5.1) 

where 

Fj(xi, ...,Xi) = /(min|xi - x,-|) (5.2) 

j<i 

with a suitable function / constructed from the zero energy scattering solution. Although the trial 
function (|5.ip is not symmetric with respect to permutations of the variables it nevertheless gives an 
upper bound to the energy of the Bose gas, because the lowest energy of the Hamiltonian is obtained for 
a symmetric wave function, cf., e.g., [25j . 

The lower bound [5] has the following main ingredients: 

1) A lemma [23] that allows to replace a possibly 'hard' interaction potential by a 'soft', nearest neighbor 
potential at the cost of increasing the range and sacrificing some kinetic energy. 

2) Division of the box into smaller boxes with Neumann boundary conditions. The size of the small boxes 
stays fixed in the thermodynamic limit. 

3) First order perturbation theory applied to the 'soft', nearest neighbor potential and an estimate of the 
error by Temple's inequality. 
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The original proof of the lower bound [5] that was limited to nonnegative interaction potentials has 
recently been extended slightly to allow also potentials with a weak negative part . 

Consider next the case of an external trapping potential V [101 I24| . For an upper bound an adequate 
trial function has the form 

N 

#( Xl ,...,x JV )=^(x 1 ,... ) x w )n^ GP ( x j) (5-3) 

with F as in (|5.1|) and i/j gp a minimizer of the GP energy functional. 

For the lower bound one first uses the GP equation to replace V by a potential depending on i(j gp . 
One approximates this potential by step functions in Neumann boxes and uses the lower bound derived 
in [5] in each box. 

In the nonrotating case this technique can be easily extended to the TF limit, g — > oo but still pa 3 — ► 

®. 

5.2 Rotating trap with g and Q fixed 

For the upper bound an apparent problem is that trial functions of the form (|5.3p with nonsymmetric F 
might lead to a lower energy than symmetric ones because the Hamiltonian is not real if f2 ^ 0. As noted 
in [llj . however, functions of the form (|5 ,3|) give an upper bound to the bosonic ground state energy as 
long as F is a real function. In fact, using the GP equation for ip GP one has in this case 

(%HV>) = NE GP +AirgN J \^ GP \ 4 

J2\ViF\ 2 -^? P \ 2 F\>+Y,Vii\F\ 2 ) IH GP (5-4) 

i<j f i 

with tp GP = ip GP (xi) and v%j = u(xj — xj). Since the right hand side of <|5 is a real quadratic form it 
does not matter whether the infimum is taken over symmetric or nonsymmetric F and a trial function 
of the form (|5.1[) can be used. 

For a lower bound the lemma of [23] is still important to 'soften' the interaction potential, but the 
division into Neumann boxes fails and an entirely different technique technique is required. This long 
standing problem was solved in [T^] by using coherent states. See also [37] for a resume. 

There remain two cases, however, that are not covered by the proof in |12j that is restricted to fixed 
g and Q: 

• The TF limit, g — > oo and/or Q — > oo, in rotating, anharmonic traps. 

• The limit — > f2 c in harmonic traps, where a GP description may completely break down. 



6 The TF Limit of the Many-Body Ground State Energy in 
Anharmonic Traps 

This section is based on [S] and concerns the leading order contribution to the many-body ground state 
energy of a dilute, three-dimensional Bose gas in the TF limit g — * oo. The rotational velocity may 
at the same time also tend to infinity. The corresponding problem for the two-dimensional GP energy 
was discussed in Section 4 where also the subleading correction to the leading TF contribution could be 
evaluated exactly. In the many-body context the extraction of the leading order contribution is already a 
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non-trivial problem. A proof that the GP energy gives also the subleading contribution to the many-body 
energy for rotating gases if — *• oo has still to be completed. 
We recall from (|2.17[) the GP energy functional 

^M=f {|(iV+ A)^| 2 + (^-if7V)|V.| 2 +^| 4 }d 3 x, (6.1) 



where the dependence on the parameters g and f2 has been explicitly indicated. The confining, external 
potential V is assumed to be smooth and satisfy V(x) > for |x| = 1 and V(Ax) = A s V(x) for all A > 
with some s > 2. 
The GP energy is 

ifE = mf{£ fl G M: |* = 1}. (6.2) 

The infimum is, in fact, a minimum and we denote any normalized minimizer by V^n- The corresponding 
density is = |^°n | 2 - 

The TF energy functional, obtained from the GP functional by dropping the kinetic term, is 



JR3 I 4 



r l p + gp z \ d d x. (6.3) 



It is defined for nonnegative densities p(-), and the TF energy 



9 

is attained for the unique minimizer: 



Effi = inf {£ TF [p] : \\p\\ x = l} . (6.4) 



TF i \ _ 1 

P g ,n W - 2 g y- g ,u ■ 4 



(6.5) 



where [■]+ denotes the positive part and p^Q is the TF chemical potential determined by the normalization 

1 1 „TF II _ I 

The TF energy and density satisfy some simple scaling relations that are useful for distinguishing 
different parameter regimes. Writing x = Ax' and p(x) = A~ 3 //(x') one obtains 

S$ [p] = A" 2 J {\°+ 2 V(x>)p> - ±(A 2 0) 2 r'V + g\-\p'f} d 3 x. (6.6) 

We now distinguish two cases: 

A. Rotational effects are at most comparable to the interaction effects: We equate A s+2 with gX" 1 , i.e., 
choose 

A = g 1 '^ (6.7) 

and obtain 

3 -/(* + 3) E TF = ^TF (68) 

with 

u = g -(s-2)/2(s + S) n (69) 

and likewise 

ff 3/(s+3 VJ F o ( 5 1/(s+3) x) = (x) . (6.10) 

Hence the TF theory has one parameter, u>. The case uj — corresponds the TF functional without 
rotation. 
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B. If to — > oo the rotational term completely dominates the interaction term. In this case it is appropriate 
to use a different scaling and equate X s+2 with (A 2 S1) 2 , i.e, take 

x = n 2 /^ 2 l (6.11) 

We then obtain 

- 2V (s-2)£TF =£ ;TF (fU2) 

with 

7 EE ft-2(s+3)/( S -2) g = w -2(.+3)/(.-2) > ( g 13) 

Moreover, asu-> oo, i.e., 7 — > 0, we have 

lim = ^TF = inf _ l r 2 . x g R 3| < U ) 

while pE\ converges to a measure supported on the set M. of minima of the function W(x) = V(x) — \r 2 ■ 
As a useful fact we note that all points in M. have the same distance from the rotation axis. 
Examples: 

If V(x) = ar s + b\z\ s , then M. is a circle in the z = plane. 

If V(x) = a|a;| s + 6|y| s + c\z\ s with a 7^ b then consists of two points. 

The scaling properties of the TF theory already suggest that one should distinguish between the 
following three cases: 

• Slow rotation, oj <C 1: The effect of the rotation is negligible to leading order. 

• Rapid rotation, ui ~ 1: Rotational effects are comparable to those of the interactions. 

• Ultrarapid rotation, u ^> 1: Rotational effects dominate. 



The main result is that in all cases the TF energy, suitably scaled, is the leading term of the many-body 

tiQM 



ground state energy as N — > 00 and g — > 00, provided the gas stays dilute, i.e., condition (|2.3[) is 



fulfilled. As a measure for the density one can use A \\pjjn\\ 00, so with a ~ g/N the diluteness condition 
can be written as A^~ 2 _g 3 1| /o^f^ 1 1 00 <C 1. 

Theorem 6.1 (QM Energy Asymptotics) Assume that N~ 2 g 3 \\p£J i \\ 00 — > as N — > 00. 

(i) 7/ <? — * 00 and — > as N 00, then 

lim { 9 -V(^3 )iv - ljE QM (7V) j = ^tf (6 15) 

(ii) If g oo and cu > is fixed as N — > oo, t/ien 

Jim {.9- s/(s+3) A^^ (A) | = £™ (6.16) 

(iii) /f £1 — ► oo and cj — > oo as TV — > oo, i/ien 

^lim {fi-V(- 2 )Ar-i£;QM (iV) | = ^tf (fU7) 

The energy estimates also lead to limit theorems for the particle density n in a many-body ground 
state: 
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Theorem 6.2 (QM Density Asymptotics for u> < oo) 

Under the conditions of Theorem 4 (i) and (ii) we have 

5 3/( S+3)A ,- V QM n ( ff l/(, + 3) x ) _ p TF (x) (6<lg) 

in weak L\ sense. 

Theorem 6.3 (QM Density Asymptotics for u> — > oo) 

Under the conditions of Theorem 4 (Hi), the scaled particle density 

n«/(- 2 )jV-V^;,n (tt 2/(s ~ 2) x) (6.19) 
converges to a probability measure with support in M. 

For the proof of Theorem 16. H it is necessary to prove first a corresponding result for the GP theory: 

Theorem 6.4 (GP Energy Asymptotics for oj < oo) 

For g — * oo with u) fixed we have 

g -s/(s + 3) E GP = E TF + Q ^- (s+ 2)/2( s+ 3) _ ( g 2Q) 

This theorem is a partial extension of Theorem 14.11 to 3D and general homogeneous potentials but since 
it only estimates the order of the subleading term in the energy and does not give not its precise form, its 
proof is simpler. In fact, the lower bound is obtained simply by dropping the nonnegative kinetic energy 
term |(iV + A)^ 2 from £ GP . The upper bound is proved by using a (suitably scaled) three-dimensional 
version of the trial function (|4.15p . In particular, this trial function exhibits a lattice of vortex lines that 
partly compensate the contribution of the vector potential to the kinetic energy term. 
In the case of ultrarapid rotations, lu — > oo, we have 

Theorem 6.5 (GP Energy Asymptotics for u — > oo) 

If ft and to — > oo then 

n -2s/(s-2) E GP = E TF + Q ^,-1 + 7 2/5^ ^ ( g 21) 

where Of = fi( s + 2 )/( s - 2 ) and 7 = w - 2 (*+3)/(«-2). 

This can be proved by picking a point Xo G Ai where V — r 2 /4 is minimal and using a trial function of 
the form 

^(x) = VM X - x o) exp {iftx • (e z A x ) /2} (6.22) 

where hg is a smooth approximation of the delta function. 

We now discuss briefly the main steps in the proof of Theorem 16.11 
1. Upper bound. 

For slow to rapid rotation (uj < oo) one can use a trial function of the form (|5.3p with a real function F 
as discussed in Subsection 15.21 The bound on E GP in terms of E TF was already described; it leads also 
to a bound on the GP density in terms of the TF density. 

For ultrarapid rotations, u) — * oo, the trial function has to be modified. One takes 

N 

#(xi, . . . , xjy) = F(xi, . . . , xjy) Y[ ^(Xi) (6.23) 

i=i 
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with F as in (15.11) and 



V>(x) = VpW exp(i£(x)) (6.24) 

where p is a regularized TF density and the phase factor S corresponds to a 'giant vortex' centered at 
the origin: 

S(x) = [±r 2 n ti\ 9 (6.25) 

where [•] denotes the integer part and is the radius of the set A4q where V — fl 2 r 2 /r is minimal. 
(This set is always a subset of a cylinder.) 

2. Lower bound. 

In the case u> < oo one first uses the diamagnetic inequality 

|(V - i A)-0| 2 > |VV>| 2 . (6.26) 

Then one writes 

V-^r 2 = V-^r 2 -p + p>-2p TF +f i (6.27) 
with the TF chemical potential p = pj^ that satisfies 

p = E TF + J (p TF ) 2 . (6.28) 

We then have to bound the Hamiltonian 

i i<j 

where g' and v' are suitably scaled versions of the coupling constant and the interaction potential respec- 
tively. This is done by introducing Neumann boxes where p TF is approximately constant and using the 
basic bound of [5] for the ground state energy in each box, 

E (n, i) > 4na{n/£ 3 )(l - o(l)), (6.30) 

where n is the particle number, I the box size and o(l) — > if an/£ 3 — > and n — > oo. This leads to the 
lower bound 

-9' |(p TF ) 2 (l + o(l)) (6.31) 

for H and altogether the bound 

NE TF {l-o(l)) (6.32) 

for the energy. 

The case of ultrarapid rotations, to — > oo, is simpler than the case of finite u> since for the lower bound 
for the energy we can ignore the (positive) interaction altogether. Namely, dropping the positive kinetic 
term, we can write for any \f 

n -2s/( s -2) N -i H ^ >c + inf w ( 6 33) 

with W(x) = V(x) - r 2 /4 and 

C^= f p' TF (x) (w^-mfW^^x + n- 23 ^'- 2 ^- 1 Yl ( 6 - 34 ) 

R3 ^ K 3 / l<i<J<.V 

Since the interaction potential v is by assumption nonnegative the same holds for , and inf W = Eq F . 
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7 Remarks on Rapid Rotation in Harmonic Traps 

In a quadratic trapping potentials, V(x) = if^ sc |x| 2 , interesting effects are expected to occur when fi 
approaches the critical frequency 

n c = V2n osc (7.1) 

at which 

y eff (x) = V(x) ~ ±fiV (7.2) 

is no longer bounded below [25], see also [H] for a review. As Sl c — O — > 0+ the effective radius R of the 
condensate tends to infinity and the system becomes effectively two-dimensional. Equating the potential 
energy corresponding to V e s and the interaction energy leads to a simple estimate for the radius: 

/ \ 1/4 

\(n c -n)nl /2 J 

The expected number of vortices is ~ flR 2 so the ratio between the particle number and the number of 
vortices is 

The expectations are now: 

1 /2 

• If AT(1 - O/nj/aOc' -> oo, the many-body ground state is still well described by GP theory. 

1/2 

• If iV(l — Q,/£l c )/aQ, c stays small, the many-body ground state is highly correlated and there is no 
GP description. The wave function should be well approximated by a function of Laughlin type in 
the variables £ = x + iy: 

*(Ci, • • • ,Cn) ~ lite - OO 2 n^C-niCilV*) (7-5) 

i<j i 

corresponding to a Fractional Quantum Hall Effect (FQHE) in the lowest Landau level. 
So far, neither of these expectations has been rigorously proved. 



8 Conclusions 

For rapidly rotating, dilute Bose gases in anharmonic traps the leading energy and density asymptotics of 
the many-body ground state as the coupling parameter g ~ Na tends to infinity can be calculated from 
the simple density functional (|6.3|) for arbitrary rotational velocities. Moreover, within two-dimensional 
Gross-Pitaeveskii theory also the subleading terms in the energy, corresponding to a lattice of vortices, 
can be evaluated exactly. The following problems are still open: 

• Prove that GP theory provides the subleading term of the many-body ground state energy for 
rapidly rotating, anharmonic traps in the TF limit. 

• For harmonic traps, prove that GP theory still applies in the limit £1 — > f2 c , provided N(l — 
Q/Q, c )/aQ^ 2 — > co. 

1/2 

• Understand rigorously the transition into the FQHE regime if N(l — Vl/Vlc) / a£l c ' remains finite. 
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